Persistent Currents in Charge-density Wave Systems 
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The inductive exchange of carriers between closed Fermi surface sections subject to Landau quanti- 
zation and open Fermi surface sections subject to charge-density wave (or spin-density wave) forma- 
tion is shown to give rise to persistent currents. This mechanism may explain recent experimental 
data in certain organic conductors in high magnetic fields. 



It is well known that a gradient in the orbital mag- 
netization within a bulk metal corresponds to a net ef- 
fective electrical current density j = V x M j^]. Such 
currents can be true macroscopic flows of charged quasi- 
particles, created in response to changes of magnetic field 
{i.e. dH/dt) and, in the case of superconductors, applied 
magnetic field H. They are considered to be "persis- 
tent" if they exist in equilibrium, and truly persistent 
currents are understood to exist only in superconduc- 
tors [||. In the type I superconducting phase, the Meiss- 
ner effect screens all electromagnetic fields from within 
the bulk 1^, while in the type II superconducting phase, 
net currents that permeate the bulk occur transverse to 
a vortex density gradient held in place by vortex pin- 
ning Long-lived currents |^ also occur in quantum 
Hall systems orthogonal to a chemical potential gradient 
V/x 1^. The equilibrium Hall electric field is sustained 
owing only to the absence of quasiparticle scattering pro- 
cesses orthogonal to the current Q. 

In this paper, I identify a further mechanism for 
persistent currents involving the coexistence of orbital 
magnetism with a charge-density wave (CDW) or spin- 
density wave (SDW) modulation With the pinning of 
a CDW (or SDW) modulation Q , a situation can be re- 
alised whereby a differential chemical potential gradient 
2V(A^) exists between the CDW (or SDW) and addi- 
tional ungapped sections of Fermi surface present. Cur- 
rents result within the bulk if, and only if, the derivative 
9M/9(A/i) is finite within the bulk. This quantity is van- 
ishingly small in most metals at low magnetic fields 
but can become quite significant when Landau levels are 
formed at high magnetic fields 0. The development of 
the gradient 2V(A/i) relies on there being a continual rel- 
ative exchange of carriers 20(AA^)/0|H| between Fermi- 
surface sheets as the magnetic field is swept [0. Again, 
this is a negligible effect at low magnetic fields, but be- 
comes very significant when Landau levels are formed at 
high magnetic fields 0. The mechanism for persistent 
currents in CDW (or SDW) systems, described here, will 
therefore be seen to be quite different than that predicted 
by Frohhch g. 

A multiband metal consisting of separate one- 
dimensional (ID) and two-dimensional (2D) Fermi sur- 
face sections is an ideal system in which to model these 



effects. Fermi surfaces of this topology are found to oc- 
cur in several organic metals and effects consis- 
tent with the existence of persistent currents have been 
observed experimentally pT|-p^. 

I begin by considering a metal that consists of only 
a single 2D pocket, characterized by a de Haas- van 
Alphen (dHvA) frequency F and cyclotron frequency 
LJc — eB/m* , where B — B • z. When the quasipar- 
ticles move in planes perpendicular to the unit vector 
z. Landau levels with eigenvalues e = hu}c{n 4- ^) are 
formed, where n = 0, 1, 2,.... The chemical potential 
/i is pinned to the highest occupied of these in order to 
conserve the total number of 2D carriers, Af2D fl^- It 
is instructive to consider contrasting situations at half- 
integral and integral Landau level filing factors v — F/B\ 
at half-integral filling factors (when -1- ^ is an integer) 
the highest occupied Landau level is exactly half filled, 
while at integral filling factors (when v is an integer) /i is 
in the transitional region between filled and empty levels. 
In the low temperature limit (T — s- 0), the corresponding 
susceptibility components Xzz — dM/dH are 
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respectively, where pxy is the off-diagonal component of 
the resistivity tensor, M — M-z and H — H-z. Note that 
the susceptibility at integral filling factors is extremely 
sensitive to Landau level broadening that results from a 
finite relaxation time r Q . It follows that M and p are 
related via the differential 
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where both e and p^y are negative for electrons and pos- 
itive for holes ]T5[ |. 

When ID carriers are added to the system, they pri- 
marily function as a charge reservoir . At half-integral 
filling factors, carriers are transferred from the ID sheets 
to the 2D pocket as H is increased, causing the magne- 
tization to increase more quickly than in Equation (Q). 
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At integral filling factors, the inverse situation becomes 
true. The presence of ID states in the Landau gaps pre- 
vents M and /x from jumping abruptly between Landau 
levels. The susceptibilities therefore become 
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where giu and (72D represent the density of states of the 
ID sheets and 2D pocket respectively 0. 

The process by which the ID sheets function as a 
charge reservoir becomes perturbed upon opening of a 
gap. For the purpose of modeling these perturbations, 
it is convenient to consider a simplified system in which 
the Fermi surface nesting gaps the entire ID density of 
states and the effects of Zeeman splitting of the ID bands 
are negligible. The latter is true in SDW systems j|] and 
also in CDW phases in which the spin up and spin down 
bands nest independently 17|. By far the most im- 
portant consideration is that the CDW (or SDW) phase 
remains stable over the field range where the dHvA oscil- 
lations occur. Persistent currents are anticipated only in 
the case of an incommensurate CDW (or SDW) in which 
the nesting vector Q can adjust itself in order to min- 
imise the total free energy of the system. They are not 
expected in the case of a commensurate CDW phase that 
is pinned to the crystal lattice . Following the solution 
of the gap equation obtained by Maki and Tsuneto , 
adapted to include the effects of a shift in the chemi- 
cal potential (labelled here as A/i) in Reference the 
change in free energy in forming such a CDW (or SDW) 
groundstate can be written as 
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The first two terms of Equation (Q) are minimised when 
the zero-temperature gap 2Ao that opens on the ID 
sheets is centred about /i, so that an offset A/i of the 
gap relative to fi costs energy. The third term, i^2D(A/^), 
accounts for an additional (either positive or negative) 
change in free energy of the 2D Landau level levels in- 
curred by this offset. The offset translates directly into 
an increase in the number of states AiViD = (/loA/i ac- 
commodated by the ID sheets, balanced by an equal and 
opposite reduction in the number of states contained in 
the 2D pocket. If Q is the component of the optimum 
nesting vector Q parallel to the mean Fermi velocity Vp 
of the ID sheets, A/i also results in a shift in this Q given 
by AQ/A/i — Qgiri/NiY), where A'^id is the equilibrium 
number of ID states. 

Should Q be uniform throughout the sample as a result 
of cooling it in the presence of a constant magnetic field. 



a finite value of A// may already be required in order 
to minimise Equation (|^) owing to the F2d(A/x) term. 
Such an effect has already been considered to explain the 
suppression of oscillations in /i observed experimentally 
in certain organic metals [Q. I will not consider this 
effect in detail for two reasons: first, it is not expected to 
operate at half-integral and integral filling factors where 
A* = -p2D(~A/i) = dF2j^{—Afj.)/dfi = 0: second, I reveal 
below that it is unlikely that a uniform A/i is maintained 
throughout the entire sample. 

If the sample is field-cooled at /i = A/i = 0, the start- 
ing condition of the sample corresponds to the point la- 
beled '0' in Fig. 0. If the CDW (or SDW) is pinned, 
this pinning initially prevents Q + AQ from changing in 
response to a change in the applied magnetic field AH. 
The optimum value of Q does, however, change, leading 
to relative difference AQ between the actual and optimal 
values, and finally to a shift in A/i. Since the density of 
states at /i in the gap on the ID sheets is zero, the sample 
begins by responding like a 2D metal in which there are 
no ID states. From this point on, the magnetization be- 
comes hysteretic. Consequently, the irreversible change 
in the magnetization is characterised by an irreversible 
susceptibility given by the difference between Equations 
(|l|) and (||); hence 
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As expected for magnetic hystersis, Xirr is always nega- 
tive. 

For a sample without boundary conditions, the mag- 
netization would continue to change with the suscepti- 
bility given by Equation (^), giving rise to the dotted 
line in Figure |^. This would occur until A/i reaches 
the thermodynamic limit, A/iiim = ± Ao/\/2, where- 
upon Equation (^ is no longer negative. At this point, 
the magnitude of the irreversible magnetization would 
saturate at a maximum value, Miir iim = ± A/iiim x 
9M/9/i ~ ± Aoz/^/2pxy, without any currents being 
induced within the bulk of the sample. Just as in super- 
conductors ^ and quantum Hall systems 0|, however, 
a net magnetization cannot exist within the bulk with- 
out incurring a current j immediately inside its surface. 
Here, the current is carried by 2D carriers, and wherever 
these exist, 



J = -V(A/i) X — , 
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where — V(A/i) is the chemical potential gradient expe- 
rienced by the 2D pocket. In order to maintain the equi- 
librium condition V/i = 0, an equal and opposite charge 
polarization field V(A/i) must also form on the nested ID 
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Fermi surface sfieets. Since CDWs (and SDWs) are elec- 
trostatically limited, rather than increasing indefinitely, 
V(AyLt) must eventually be limited by the characteristic 
threshold electric field 



eEt, 
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required for its depinning. This equates to a critical cur- 
rent density 
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for the carriers in the 2D pocket within the sample. 

From Equation (||) , one can see that the electrodynam- 
ics has parallels with the quantum Hall effect Q, with 
the exception that there exists zero net electric field in 
equilibrium: the quantity Et instead parameterizes the 
charge polarization field required to slide the CDW (or 
SDW) , that ultimately leads to dissipation. The absence 
of a net electric field implies that the effect, described 
here, is not limited only to integral Landau level filling 
factors. The inductive behaviour of the sample is, in 
contrast, much more like that of a type II superconduc- 
tor, but with an irreversible susceptibility Xin that can 
depart significantly from the ideal diamagnetic value —1. 
As AH increases, the surface region in which Et is polar- 
ized propogates further into the sample. The inductive 
response of the sample to a change in applied magnetic 
field AH is therefore indistinguishable from the critical 
state model developed by Bean (2l|, provided Et is fi- 
nite for all directions within the planes. This is certainly 
the case in CDW (or SDW) systems for which Q occurs 
at an oblique angle with respect to the lattice vectors, 
where Fermi surface nesting results in charge (or spin) 
modulation that is two-dimensional. 

If one considers a cylindrical sample of radius r, with z 
aligned along its axis, following Bean the irreversible 
magnetization becomes 



Mirr = -Xii 



AH- 



{AHf {AHf 



H* 



3H 



*2 



(9) 



This region is identified by point '1' in Fig. |^. In this 
model, H* is the characteristic coercion field while Xm-H* 
is the maximum magnetic field that can be screened 
within the sample. When AH = H* , corresponding to 
point '2' in Fig. 0, the irreversible magnetization reaches 
its saturation limit for a cylindrical sample. 
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When the direction in which the field is changed is re- 
versed at point '3' in Fig. Q the initial susceptibility 
is, once again, given by Equation (^. A similar cubic 
law should apply to that in Equation (|^) in the reversal 
state, but with double the interval in AH being required 



to reverse the polarity of V(A/i). The magnetization sat- 
urates again at point '4,' and a full hysteresis loop results 
on further cycling of AH. 

I now compare the predictions of this model with re- 
cent measurements made on samples of the organic con- 
ductor a-(BEDT-TTF)2KHg(SCN)4 at high magnetic 
fields. This material has a Fermi surface which nearly 
matches that of the model [|lO|. At low temperatures, 
T < 8 K, and low magnetic fields, figH < 23 T, a- 
(BEDT-TTF)2KHg(SCN)4 is believed to possess a com- 
mensurate CDW phase, transforming to one that is in- 
commensurate at high magnetic fields hqH ^ 23 T 
[|ll|-||,|7|jl|j2|-||l. It is within the high magnetic field 
phase where effects consistent with persistent currents 
have been reported |11-14|. 

I begin my comparison by considering the irreversible 
susceptibility. On computing Xirr, using Equation ^ 
and the established parameters {F « 670 T, m* « 2 me, 
N2B ~ 1.6 X 10^6 m-3, « g2D and T « 2 ps || for a- 
(BEDT-TTF)2KHg(SCN)4), the model estimates (xmod) 
listed in Table | are found to agree rather well with the 
experimental values (xoxp) [0 both at half-integral and 
integral filling factors. 

I now estimate the threshold electric field from the 
saturation magnetization Miir^max (listed in Table Q as 
Afmax)- This is found to be marginally higher at half- 
integral filling factors than at integral filling factors in 
Reference but becomes comparable at the lowest 

temperatures |jl^. Making the cylindrical sample ap- 
proximation, with r ~ 0.5 mm, one can estimate the 
mean scalar critical current jc using Equation ( p^ ) and 
the corresponding threshold electric field Et using Equa- 
tion (||). I believe the estimates in Table | to be physi- 
cally realistic because the total potential difference across 
the sample Vt = rEt does not exceed the maximum 
thermodynamic potential V2eAo ~ 0.7 meV (where 
Ao - 0.5 meV U) that can be sustained without de- 
stroying the CDW groundstate. 

Next, I consider the coercion fields required to reverse 
the polarity of the current. On applying the simple ex- 
pression for H* in Equation (^|) , the values estimated in 
Table l] compare favourably with experimental observa- 
tions [l3. The behaviour of the irreversible magnetiza- 
tion is also found to be very well described by a cubic law 
of the form given by Equation (H). Note that twice the 
interval in field 2H* is required to reverse the polarity of 
the currents throughout the entire sample. Meanwhile, 
the maximum magnetic field that can be screened by the 
sample, XmH*, is a small fraction of a millitesla. 

Finally, I consider an odd situation that occurs at 
the minima and maxima of the dHvA oscillations where 
the orbital susceptibility x = dM/dH vanishes. Con- 
sequently, there will be discrete fields at which Xm ~^ 
causing the irreversible magnetization never to reach sat- 
uration within the quarter period of the dHvA oscilla- 
tions AH = hqH^ /AF. Because there exists a chemical 
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potential gradient V(A/x) across the sample, however, 
the situation where Mi„ = can only exist in a narrow 
region of annular topology within the interior of the sam- 
ple. The effect should still be observable, however, as 
drop in the volume-averaged saturation magnetization. 
This could provide an explanation for the drop in the 
saturation magnetization observed at intermediate filling 
factors at the lowest temperatures in Reference 

In conclusion, I have identified a mechanism for persis- 
tent currents in metals in which a CDW (or SDW) phase 
coexists with a strongly field-dependent orbital magneti- 
zation. The magnetic field drives an exchange of carri- 
ers between the ID and 2D Fermi surface sheets, which 
has the effect of stretching or compressing the CDW (or 
SDW) like a concertina. Persistent currents result when 
pinning of the CDW (or SDW) resists the concertina- 
like motions. The charge polarization field of the pinned 
CDW negates the Hall electric field that would normally 
give rise to the Hall effect Q]. However, because the 
CDW (or SDW) exists in a critical state on the brink of 
sliding, the magnetic properties resemble those of a type 
H superconductor |21|. Like a superconductor, currents 
exist only because of impurity pinning. The mechanism 
for persistent currents, described here, however, is nei- 
ther the quantum Hall effect or superconductivity, but 
one that is entirely different. I propose this effect as an 
explanation for persistent currents observed in a-(BEDT- 
TTF)2KHg(SCN)4 
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FIG. 1. Magnetic hysteretic behaviour anticipated by the 
persistent current model for a cylindrical sample, as described 
in the text. 
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TABLE I. A comparison of the irreversible susceptibility 
of a-(BEDT-TTF)2KHg(SCN)4 expected for an ideal sample 
according to the model, Xmod, with that obtained experimen- 
tally, Xcxp [12]. Also listed are the maximum irreversible mag- 
netization, Mmax, obtained experimentally together with the 
critical current density, jc, CDW depinning threshold electric 
field, Et, total potential difference, vt, coercion field, fioH* , 
and screening field, /xoXexp-ff*, calculated according to the 
model. 
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Figure 1 of Harrison et a! 



